Casimir energy in finite expanding universes: a prelude 
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We present the Green's functions that are the solutions of the massive Klein-Gordon equation for 
a scalar field with non-minimal coupling to gravitation for several static and expanding cosmological 
models. An important feature of such functions is that they can be used to study the appearence 
of a Casimir energy in models of the universe with compact spatial sections. 
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I. INTRODUCTION 

In the description of the universe based on General 
Relativity, the idea of curved manifolds is essential. Ein- 
stein's equations are mathematical relations between the 
curvature of spacetime and the matter and energy, which 
act as a source of the curvature and can be seen as a 
first step in the study of the structure of the physically 
possible models of the universe. Given the curvature of a 
spacetime, it is important to characterize its global prop- 
erties, such as its topology. Since General Relativity pro- 
vides no clues about these properties, it is valid to include 
in the list of possible models of the universe, those in 
which the space does not have the simplest topology. 

A space with non-trivial topology may be compact 
(i.e., finite), regardless of its curvature. This means that 
the universe may present certain peculiarities due to its 
finite size, such as periodic spatial patterns. The observa- 
tion of any of these peculiarities may indicate the global 
topology of the universe, so that it is very important to 
know what kind of mark a determined topology may im- 
print on the visible universe. Cosmic topology, which 
deals with this issue, has been the subject of a grow- 
ing interest, reflected in the publication of recent reviews 
[1, 2] and introductory texts [3, 4, 5]. 

An interesting physical effect due to the finiteness of 
the space is the appearance of a Casimir energy. The sim- 
plest case where this energy, measurable experimentally, 
appears is between two parallel reflecting planes in the 
Minkowski vacuum. However, this case can be general- 
ized. The Casimir energy can be seen to be created by the 
imposition of periodic boundary conditions on a generic 
field propagating in the vacuum and can be studied theo- 
retically with a scalar field, by means of Green's functions 
G(x,x') for the Klein-Gordon equation [6]. The proce- 
dure begins with the classical energy-momentum tensor 
associated with the field, obtained from the Green's func- 



tion through a second order differential operator 



T^iv {x, x') — -D^jf G (x, x ) . 



(1) 



From this, we obtain the mean value of the field in the 
vacuum, such that 



<0|T„„|0) - Urn T„ v (x, x') . 



(2) 



In general, when speaking of the Casimir energy, it is to 
the finite part of the 00 component of this mean value 
that one refers (see, e.g., [7, 8]). 

As seen in some recent works [7, 8], the same kind of 
study can be done in models of the universe with non- 
trivial topology. However, these studies have been lim- 
ited, until now, to static models of the universe. Never- 
theless, there is not, a priori, any reason not to include 
models with expansion. What is important to notice 
is that a fundamental step in this process is to write, 
for each kind of expansion of the universe, the specific 
Green's functions for the Klein-Gordon equation. 

The purpose of this work is to present the Green's func- 
tions found as solutions of the generalized Klein-Gordon 
equation for a scalar field of mass m with non-minimal 
coupling £ to the gravitational field for static and ex- 
panding cosmological models. In the expanding cases, 
the functions presented are valid for a group of cosmo- 
logical models, parametrized by the curvature scalar k 
and by an index n. Such a group of models includes, for 
example, the de Sitter solution. 

The Green's functions for universes with non-trivial 
topology are obtained from the general solutions by not- 
ing that the Green's function for a generic finite space, 
M = M/T, is 



G M (x,x') = J2 G m 



x,T(x') 



(3) 
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where T are the elements of the group which defines M. , 
each one represented as an operator that identifies a point 
P on the finite space M. and its image P' onto the infinite 
covering space M, i.e., P 1 = T (P) = P, with P' e M 
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and P £ M. . For 3-dimensional spaces of constant curva- 
ture the covering space M can be [1] 

• the spherical space S 3 , of positive curvature; 

• the fiat space E 3 , of null curvature; 

• or the hyperbolical space H 3 , of negative curvature. 

The paper is organized as follows. The Klein-Gordon 
equation for a scalar field is first written in a general- 
ized form valid for the curved spacetimes described in 
the Friedmann-Lcmaitre-Robertson- Walker line element 
of standard cosmology. Some general Green's functions 
are then obtained for specific cases of static and ex- 
panding spacetimes. Finally, some comments on the re- 
sults obtained are given. General solutions of the gen- 
eralized Klein-Gordon equation are presented in the ap- 
pendix. Throughout this work, natural units are used 
(c = G = h=l). 



for the evolution of the scale factor a, with p the energy 
density of the matter, and 



1 d 2 a k 
a 6 dij z a 1 



(9) 



where p is the pressure of the matter. For the FLRW 
metric, the Ricci scalar is 



R = 



I (Pa 

a drj 2 



+ k 



and 



□ = 



d 2 d 2 

drj 2 drj 



Therefore, the field equation becomes 



1 



d 2 <p 

drj 2 



+ 2D^ - VV + (£R + m 2 ) a 2 if 



(10) 



(11) 



0. (12) 



II. INITIAL CONSIDERATIONS 



III. GREEN'S FUNCTIONS 



For a classical scalar field with a Lagrangian density 



one obtains the field equation 

Dip + (£R + m 2 )<p = Q 



(4) 



(5) 



where □ = [—g] <9 M (\/— gg^dv) is the d'Alcmbertian 
operator, m 2 and £ are constants, and R — g llv R ,lv is the 
Ricci scalar derived from the spacetime metric g^ v , with 
g = det [g^], 

In this work the interest is focused on metrics repre- 
senting cosmological solutions, in particular, those with a 
Fricdmann-Lemaitre-Robertson- Walker (FLRW) line el- 
ement, 



ds 2 = dt 2 -a 2 (t) 



dr 2 



1 — kr 2 
which can be written as 



ds 2 = a 2 fa) 



d v 2 ~d x 2 - 



sin 



k 



+ sin 2 9d(j) 2 ) 



I 2 + sin 2 



(6) 



(7) 

where k = 0, ±1, with t the cosmological time, r] the con- 
formal time, and dt — adrj. Assuming Einstein's equa- 
tions with a cosmological constant A, and matter a per- 
fect fluid, this line element gives the Friedmann equation 



D 2 1 (lda\ 8tt A k 

3 = ^ :x =^p+o-^ ( 8 ) 



a 2 \a drj 



A Green's function can always be constructed, using 
the eigenfunctions and eigenvalues of the general solu- 
tion of the homogeneous equation. However, this is not 
always an easy task, especially in non-flat spacetimes, 
where, for example, plane waves are replaced by general- 
ized solutions to Helmholtz equation [9] (see appendix). 
Therefore, since the Green functions studied here must 
satisfy the more general equation, 



1 



8 2 G 



drj 2 



2D 



dG 



V 2 G 



+ (£R + m 2 ) G = 



S(x- x') 



(13) 

we search for solutions to Eq. (13) in its homogeneous 
form, keeping in mind that such solutions must be sym- 
metric in x and x' . In the appendix, some general solu- 
tions of the homogeneous field Eq. (12) are presented. 

Assuming isotropy of the universe, the use of the coor- 
dinates (x, @, 4>) allows for a choice of orientation where 
the two points x and x' are radially aligned, i.e., 6 = 8' , 
4> = </>', and x' = 0, so that for the FLRW line element, 



d 2 d 

V 2 = — + 2VkcotVk X jr> 
ox ox 



(14) 



where x is the spatial distance between the points x and 
x'. Thus, the homogeneous equation to be solved is 

d 2 G n dG d 2 G 2Vk dG 



drj 2 



d-q dx 2 tan Vk~x d X 



-2D— — - + (£R + to 2 ) a 2 G = 0. 

(15) 

Several different Green's functions obey this equation 
since there are several boundary conditions that can be 
defined. For example, according to different boundary 
conditions, we have the causal Green function (or Feyn- 
man propagator) Gf{x,x'), the Wightman's functions 
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G + (x, x') and G~ (x, x') for positive and negative fre- 
quency, respectively, or the elementary Hadamard func- 
tion G^ (x, x 1 ). In quantum field theory, these different 
functions are associated with mean values of some field 
functions [6]. 



A. Static universes 

For flat spaces, the trivial solution of Fricdmann's 
equation is Minkowski space, which is a vacuum static 
solution for A = 0. In this simple case, a = 1, R = 0, and 



d 2 G d 2 G 
drj 1 d\ 2 

Using the invariant 



2 dG 2 ^ 
_ —+m 2 G 
X ox 



0. 



VArf 2 



X 



where Arj = 77 — 77', Eq. (16) becomes 

d 2 G | 3dG | m 2 G _ 
ds 2 s ds 

with the general solution 



1 



G (x, x') = - Ci h[ 1} (ms) + c 2 H^> (ms) 



(2) 



(16) 



(17) 



(18) 



(19) 



where and are Hankel's functions of order v 

of the first and second kinds, respectively [10]. Another 
invariant, the world function a [11, 12, 13], which in the 
Minkowski space is given as a = s 2 /2, could have been 
used to give a similar general solution. Values for the coe- 
ficients ci and C2 are set in accordance with the particular 
boundary conditions used in different Green's functions. 
For example, in the case where k = m — 0, Eq. (13) 
is a wave equation with the massless positive frequency 
Wightman function, 



D+ (x,x') = - 



4ir 2 s 2 ' 



(20) 



Another type of Green's function, the massless 
Hadamard elementary function, 



(x,x') 



1 



(21) 



is also a solution to Eq. (13) for this case. Since near the 
origin [10], 



H P(x)*- — , HY>{x)~—, 



2i 



(2) 



2i 



nx 



nx 



we have 



G + (x,x') = ^ 2) (ms). 



(22) 



(23) 



This solution can be used as a limit that any other posi- 
tive frequency solution for curved spaces must reach. 



For non-flat spaces (k = ±1), Friedmann's equation 
with a non-zero cosmological constant has the vacuum 
(p = 0) static solution, 



a=(3fc) 1/2 A- 1 / 2 , 



(24) 



with kA 1 > 0. In this case, we have, from Eqs. (10) 
and (12), 



d 2 G d 2 G r- rr dG 

JT2 ~ 7TT ~ 2 vfccot Vkxjr- 
dr/ 2 d\ 2 ox 



6C 



3m z 



where k = ±1. If 

G(x,x') 



Vk~X 



■A(x,x') , 



sin-\/fcx 

Eq. (25) becomes similar to (18), and thus 
Vkx 1 



kG = 0, 
(25) 

(26) 



G (x, x) = 



and 



sin Vk~x s 
G+ (x,x') 



Cl H^> (ps) + c 2 HY> (us) 



(2) 



(27) 



* ^ X ^f) M , (28) 



87T s i n y/kx s 

where s is given by Eq. (17), and n 2 = 3m 2 A -1 + 
(6£ — 1) k. Therefore, when 6£ = 1, m = 0, and x — 0, 
we again obtain Eq. (23), as expected. 

The spherical case, k = 1, deserves special considera- 
tion since the universe represented is finite spatially, so 
that there are several equivalent spatial paths between 
the points x and 2/, all of which must be represented 
by the Green's function. Therefore, using the internal 
symmetry of the space § 3 in Eq. (28), we have 



G+ {x,x') 



8tt 



v 2 - xl 



smx 



, (29) 



where sin^n = sin (x + 2nn) — sinx was used. For any 
other spatially finite universe one uses an equivalent pro- 
cedure [7, 8]. 



B. Expanding universes 

When A = 0, a group of solutions of the Friedmann 
equation (8) can be found by means of the ansatz 



p = Ca- n , 



(30) 



where C and n are constants. This general group of so- 
lutions is then 



a (r)) = a 2 -™ /c 2 -- esc 2 -" \fk (J^ - l^J 77 
where a 2 = 8ttC/3 and n ^ 2. Thus we can write 
L> = Vfccot Vk (| -l) 



(31) 



(32) 
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n 


a[t) 


Comments 





(2a) [e + ke J 


exponential expansion 


1 


2- 2 a 2 t 2 + ka-' 2 


quadratic expansion 


2 


\JcP- — kt 


linear expansion 6 


3 


2- 2 ' 3 (3a) 2/3 i 2 / 3 


matter" 


4 


V2at - kt 2 


radiation 



"Exponential expansion is characteristic of the de Sitter's solution, 
a vacuum solution with a cosmological constant. 

^Linear expansion is characteristic of Milne's universe, where k = 
— 1 and p = 0. 

c Only the solution with k = 0, known as Einstein-de Sitter, is 
shown here. The more general solution can only be written as a 
function of the conformal time r). 

TABLE I: Summary of some of the solutions (not all phys- 
ically relevant) of Friedmann's equation with A = for the 
scale factor a (£), from the ansatz p = 3q 2 (87t) _1 a~ n . 



and 

Ra 2 = (4 - n) 3k esc 2 Vk Q - l) rj . (33) 

Using the ansatz given by Eq. (30) we can solve Fried- 
mann's equation for a and then substitute the expression 
found into Eq. (9), to obtain the linear equation of state 
3p = (n — 3) p. 

For the case where n — 2, the solution is 



a (rj) = cxp 



(a 2 -k) 1/2 77 



with 



D=(a 2 -k) 1/2 ; Ra 2 = 6a 2 . 



(34) 



(35) 



Table I shows, for various integral values of the index 
n, the expression for the scale factor a as a function of 
the cosmological time t. 

Although Eqs. (31) and (34) give a group of solutions 
which are valid for any value of the parameter n, not 
all values of n give easily solvable field equations. Here 
we consider only the cases where n — 0, 2 and 4. For 
n = 2 and n = 4, only the massless (m = 0) equations 
are considered. 



1. Solution with n — 

Since in this case p is a constant, one has then a situa- 
tion completely equivalent to the case of a vacuum with 
the presence of a cosmological constant A = 3er 2 . There- 
fore, the solution with n = is equivalent to the de Sitter 
solution, the curved spacetime most studied by quantum 
field theorists [6], and which is to cosmology what a hy- 
drogen atom is to atomic physics [14] . 

The de Sitter solution is, in fact, a group of solutions, 
each one with a different k, but all related to portions 



an hyperboloid embedded in a 5-dimensional Minkowski 
space of metric rj^ u = diag [1,-1,-1,— 1,-1]. Different 
sections of this surface represent the various de Sitter 
spacetimes, each one with a distinct curvature. There 
also exists a static representation of de Sitter spacetime, 
which will not be treated here [6, 15, 16, 17, 18]. It 
is, however, worthwhile to note that despite much in- 
vestigation of the de Sitter spacetime, it is rare to find 
descriptions of the hyperbolic de Sitter solution, which is 
not even cited in some standard texts [6, 18]), perhaps 
due to an 'aesthetic' prejudice: "the open model has the 
unnatractivc feature that in horosphcrical coordinates it 
leads to an uneasy distinction of the third spatial com- 
ponent" [19]. 

For the flat case, we have 



ds 2 = a- 2 V - 2 [air, 2 - d X 2 - X 2 (d9 2 + sin 2 Odtf)] , 



with the corresponding parametrization 



(36) 



(37) 



x° = 2- 1 7 ? - (a- 2 + x 2 ) (2?7) 

x 1 = -2-^ - (a- 2 - x 2 ) (2J])- 1 

x 2 = — (arj) 1 xsinf)cos0 

x 3 = — (arj) ~ 1 x sm sm </> 

x 4 = — (ar/) 1 x cos 9 



where < (f> <2ir, < 6 < ir, 0<r<oo, -co < r\ < 0. 
In this case, it is easy to see that the parametrization of 
the hyperboloid is incomplete since x° + x 1 > 0. 
In the spherical case, the line element is 



ds 2 = a 2 esc 2 rj [drj 2 — dx 2 — sin 2 x {d0 2 + sin 2 



with the relations 



x° = a 1 cot rj 



x = oT L esc 77 cos x 
x 2 = aT 1 esc 77 sin x sm & cos < 
x 3 = aT 1 esc 77 sin x sm # sin <! 
x 4 = aT 1 esc 77 sin x cos 



where O<0<2tt, 0<6 <n, 0<x<tt, < 77 < tt. 
Finally, for k = —1, we have 



(38) 



(39) 



ds 2 -- 
with 



csch 77 



[dry 2 - d X 2 - sinh 2 X (d0 2 + sin 2 8d<j) 2 



of the 4-dimcnsional surface 77^^^^ 



-a , which is 



(40) 



(41) 



where < (f> < 2n, < 6 < n, 0<x<oo, —00 < 77 < 00. 

The invariant quantity to be considered in this case 
is not simply the interval s, as defined in Eq. (17), but 



' x° 


= or 


1 csch?7COshx 


x 1 


= or 


1 coth 77 


x 2 


= or 


^xsch 77 sinh x sin 9 cos 


X 3 


= or 


1 csch 77 sinh x sin 9 sin 


U 4 


= or 


1 csch?7 sinh x cos 6 
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rather the geodesic interval a(x,x') in the hyperboloid 
[20, 21], a geometric quantity which is another example 
of a world function [11, 12, 13], given in this case as 
[22, 23, 24, 25] 



a(x,x') = a arccosh \p (x, x')\ , 



(42) 



where p (x, x') = —a 2 i]^ u x^ (x') v is an auxiliary quantity. 
Using Eqs. (37), (39), and (41) and the radial alignment 
between the points x and x', we can write 

' 1+(A ?? 2 - X 2 )(2 W ')" 1 
p(x,x')=l 1 — esc rj esc r/' (cos A?7 — cos x) > 
1 + csch rj csch r( (cosh At] — cosh x) 

' (43) 

for fc = 0, 1, and —1, respectively, so that, for any value 
of k, we obtain 



(P 2 - 1) 



d 2 G 



dp 2 

with the solution 



dG_ 

dp 



(12£, + m 2 a- 2 ) G = 0, (44) 



G{x,x') = Cl F Q - i/, | + i/; 2; 

^ /3 3 „ 1 +p 



(45) 



where ^ 2 = 2~ 2 3 2 - (12£ + m 2 a~ 2 ). An analysis of the 
conditions that must be imposed on this solution to ob- 
tain the values of C\ and c 2 can be found in the literature 
[25] and will not be reproduced here. 

We note that Eq. (44) has also solutions in terms of 
the associated Legendre functions [22, 23], 



G(x,x') = (p 2 - 1) 1/2 \ Cl P\ +iy (p) + c 2 Q\ k+v (p) 



or 



G(x,x') 



dp 



ciP i +v (p)+<£Q i+„(p) 



valid since p assumes only real values. 

In the specific case of a massless field (m 
conformal coupling (£ = 1/6), we get 



CiP + c 2 

p 2 -l ' 



(46) 
, (47) 
0) with 
(48) 



which gives 
D+ {x,x') = - 



8tt 2 (p - 1) 



167T . 2 aa 

— -sinh — 

a z 2 



(49) 

a result formally similar to the one obtained for an uni- 
formly accelerated observer in Minkowski space (see also 
Eq. (20)) [6]. 



2. Solution with n — 2 

For n = 2 we have the line element 
dr 2 



ds 2 = e 2r >^^ 



drf 



1 — kr 2 



- r 2 (d6 2 + sin 2 8d<j) 2 



(50) 

which, for k — — 1 and a = 0, can be transformed in 
Minkowski space by the change of coordinates 



r = e r >r, t = e v y / l + r 2 . 



(51) 



This peculiar parametrization of the flat Minkowski space 
is known as the Milne universe [26, 27]. 

In the general case where n = 2, the homogeneous Eq. 
(15) becomes simpler in the flat massless case, i.e., with 
k = m = 0, and we have 



(52) 



(53) 



d 2 G n dG d 2 G 2dG a& 2 ^ n 

— + 2a— — -^~2 -j. 1- 6?a 2 G = 0. 

dr] 1 dr] d\ z X OX 

This damped wave equation [28], has the solution 

G (x, x') = e~ aAv Go {x,x'), 
where Go (x,x') obeys the equation 

9 2 Go d 2 G 2 <9G , . 2r < n iz.a\ 

-W~l^-xlh +m - 2)aGo = > (54) 

which is Eq. (16), for Minkowski space. Therefore, 



G(x,x') = e - aAr >- 
s 



aHP ((3s) + c 2 HY> (/3s) , (55) 



(2) 



where f3 2 = (6^ - 2) a 2 and s is given by Eq. (17). The 
non-flat massless cases can be obtained from this solution 
in a similar way to that used for static non-flat cases. 



3. Solution with n = 4 

Although for n = 4, we have Ra 2 = from Eq. (33), 
it is not easy to find a Green's function in this case. 
However, it is interesting to see that in the homogeneous 
massless (ra = 0) equation there exists a symmetry be- 
tween the coordinates 77 and x, i.e., 

d 2 G n rr rr dG d 2 G n r- r- dG „.„, 
dr] 1 dr] dx d\ 



IV. FINAL REMARKS 

We presented the Green's functions for a number of 
interesting cosmological solutions. Not all plausible pos- 
sibilities for scalar fields have been covered and there re- 
mains a number of physically possible cosmological solu- 
tions, for example, the expanding solutions with n = 2 
or n = 4, for which there is not yet a complete, easy 
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available, analytical treatment. Even where many stud- 
ies have been made, there are lapses, such as the almost 
complete absence of references to the hyperbolic de Sit- 
ter solution. It is the intention of this article to stimulate 
more work in this field. 

In the context of cosmic topology, the task of writing 
Green's functions is a preparatory step, and can be done 
analytically. The bulk of the work, the use of these func- 
tions in cosmic topology, is mainly numerical and, for 
this reason, is not presented here. 

The importance of studies linked to the field of cos- 
mic topology must not be underestimated. Despite some 
recent claims about the beginning of a new era of 'preci- 
sion cosmology', we are still unable to determine whether 
the universe is finite or not. Moreover, with the ob- 
servations being made currently, we cannot tell for sure 
whether the universe has a really null curvature or just 
a very small one, which would be a remnant from a pre- 
inflation era. In this context, finding topological proper- 
ties of the universe is very important. In the theoretical 
search for these properties, the study of a scalar cosmo- 
logical Casimir effect in finite spaces, while not providing 
the final answer about the shape of the universe, would 
certainly supply important information. 
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APPENDIX: SOLUTIONS FOR THE FIELD 
EQUATION 

A Green's function can always be written using the 
eigenfunctions and eigenvalues of the homogeneous equa- 
tion it obeys [28, 29]. However, to obtain Green's func- 
tions by this process for the finite spaces of non-trivial 
topology is not an easy task, since the eigenvalues in gen- 
eral are not easy to find. Therefore, it is more convenient 
to write a Green's function for infinite space and force 
it to satisfy the symmetries of the finite space numer- 
ically. Nevertheless, even the functions for the infinite 
space are not easy to obtain in this fashion because the 
eigenfunctions are, in general, very complicated. In this 
appendix, eigenfunctions for some specific cosmological 
solutions are presented. 

The field equation (12) has general solutions of the 
kind 



where u and v are functions satisfying the equations 
d^ u du 

d^ + 2D d^ + ^ + ^ R + m2 ) fl2 ] u = °' (A - 2) 

and 

V 2 v = -~£ 2 v, (A.3) 
the Helmholtz equation. 

1. Spatial modes 

To obtain information about the size and shape of the 
universe, we must study the spatial modes of the general 
solutions of the field equation. A general form for the 
solutions of the Helmholtz equation, without taking into 
account the finiteness of the universe, is [6] 

fexpr^V, £ 2 = r ]ij £ i P (fc = 0) 

k ; \ n£> ( X ) Yf (9, J>), £ = (£, J, M) (k = ±1) ' 

(A.4) 

where Yf 1 are the usual spherical harmonics (with limi- 
tations on the possible values of £, J, M), and 



n « (X) = c/j sinh J X 



d cosh x 



l+J 



cos^x, (A.5) 



with IL^ (x) obtained by replacing I with — i£ and x 
with — %x [6]. 

We note that the solutions given above are not the only 
possible options: for 3-D flat spaces, there are 11 coor- 
dinate systems in which the Helmholtz equation is sep- 
arable [29, 30], while for hyperbolic spaces, there are at 
least 12 such systems [31]. Such diversity is particularly 
important in universes with compact spatial sections. In 
these universes the coordinate system (x, 6, <f) is not al- 
ways a good choice since the characteristic symmetries 
of the spatial sections of such universes generally do not 
fit in these coordinates, especially in the hyperbolic case. 
An example where this occurs is in the calculation of 
the volumes and other related functions of the polyhcdra 
that serve as the fundamental cells of compact hyperbolic 
spaces and the use of the coordinate set (x, 0, (j>) makes 
the problem almost intractable [32]. 



2. Non-spatial modes 

Trying a solution of the kind 

u fa) = aPU (77) , (A.6) 
where p is a constant, in Eq. (A. 2), we obtain 



tp(r),x) =u(i])v(x) , 



(A.l) 



+P 



„ „o 1 d 2 a 
(p+l)D 2 + - — 
a drj z 



U = 0. (A.7) 
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In the case where p = — 1, we obtain , using Eq. (10), with 



d 2 U 
dr] 2 



a 2 } U = 0. (A.8) 



U = Ci sin Ary + c 2 cos A?7 . 



(A.11) 



This equation, which clearly becomes simpler for the case 
of a conformal coupling, £ = 1/6, can be solved analyti- 
cally for several cases of cosmological interest. 



a. Static universes 



For Minkowski space, a — 1, k = R = and 



U = c\ sin 



(l 2 + m 2 ) 1/2 V 



+ C2 COS 



(^ 2 +m 2 ) 1/ % 



For fc ^ 0, we have from Eqs. (10) and (A.8), 



(A.9) 



= -(V + 6e + ™ 2 !)l^-A 2 C/, (A.10) 
I 



b. Expanding universes 



When A = and the scale factor a is given by Eq. 
(31), the function u satisfies the equation 



d 2 u 2D^ U 
dr/ 2 dr\ 



Then, 



3fc£ (4 - n) u 



sin 



v^(i-i) 



(l 2 + m 2 a 2 ) u = 0. 

(A.12) 



d 2 Uk,n 
dr] 2 



+ Ik 



(2 + k(A-n) (6g-l) + 



m 2 Q! 2-nfc2- 



2 sin 2 



sm 2 



\/fc(f - 1) r? 



(A.13) 



where the expression for a given in Eq. (31) was used. 



Solutions for n — 2 



Solutions for n — 
When n = 0, Eq. (A.13) becomes 



d 2 U kfl 
dr] 2 



+ 



k + £ + ■ 



2 )k 



sin 2 \fkr] 



U kfi = 0, (A.14) 



where v 2 = 2~ 2 3 2 - (12£ + m 2 ^ 2 ) . The solutions to 
this equation are 



U ,o = Vv [ciJ* (iv) + C2N U (in)] , 
for k = and 

f/i,o = sin5»j Cl p^i +v ^_ T (-cos?7) 

+ c2 ! :i +VFTT (-cos^ 

and 



(A.15) 



(A.16) 



[/_i,o = sinh 2 77 ciP^i +v ^ r - T (cosh 77) 

+ c 20":i +VF 3T(cosh77)] (A.17) 

for k = 1 and fc = — 1, respectively. 



When n = 2, Eq. (A.8) is 



d 2 U k ,2 
drf 



+ 



£ 2 + k + m 2 e 2r '^^ + (6£ - 1) . 



which has the solution 

U k ,2 = Cl J 7 

where 



U k , 2 = 0, 
(A.18) 



( a 2 -k) 1/2 v 
(a 2 - fc) 1/2 



,(a a -fc) 1/a , 



(a 2 - fc) 1/2 



(A.19) 



7 2 = -(a 2 -fc) 1 [^ 2 + fc+(6£-l)a 2 ] . (A.20) 
This solution reduces to 

C/™2 =0 = ci sin (777 \/ a 2 — fc^j + c 2 cos (^yr]\/ a 2 — kj 

(A.21) 

when m = 0. 

Solutions for n = 4 



When n = 4, Eq. (A.13) becomes 



dr] 2 



2„,2 



k + l 2 + 



■ sm 



I7*,4 = 0, (A.22) 
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whose solution for flat space is 



and a modified Mathieu equation, 



U 0A 



1 



ciM 



(imar) 



+ c z M P{um a )-\±-^{ imOL tf) . ( A - 23 ) 



where 



M AiM (z) = / +1 / 3 e-'/2$ ^ - A + 1 2 M + 1; . 

(A.24) 

is the Whittaker M function [10]. This solution can also 
be written using the parabolic cylinder function D p (z) 
[10]. 

For k = 1 and fc = -1, Eq. (A. 22) yields a Mathieu 
equation, 



drf 



t + 1 + 



m 2 a 2 m 2 a 2 



cos2?7 



C/i,4 = 0, 
(A.25) 



d 2 /7_i 4 

— + 

drj 2 

respectively. 



m 2 a 2 m 2 a 2 



} + — ^— cosh2?y 



U- 1A = 0, 
(A.26) 



In the massless case, Eq. (A. 22) has a single solution 
for all fc's, 



£/™ 4 =° = ci sin 



2\V2 



c 2 cos 



(*• 



^1/2, 



(A.27) 
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